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t Abstract Millisecond pulsars represent an evolutionar- 
' . ily distinct group among rotation-powered pulsars. Out- 
side the radio band, the soft X-ray range (~ 0.1-10 keV) 
is most suitable for studying radiative mechanisms op- 
crating in these fascinating objects. X-ray observations 
revealed diverse properties of emission from millisecond 
pulsars. For the most of them, the bulk of radiation is of 
a thermal origin, emitted from small spots (polar caps) 
on the neutron star surface heated by relativistic parti- 
cles produced in pulsar acceleration zones. On the other 
hand, a few other very fast rotating pulsars exhibit al- 
most pure nonthermal emission generated, most proba- 
bly, in pulsar magnetospheres. There are also examples of 
nonthermal emission detected from X-ray nebulae pow- 
ered by millisecond pulsars, as well as from pulsar winds 
shocked in binary systems with millisecond pulsars as 
companions. These and other most important results ob- 
tained from X-ray observations of millisecond pulsars are 
reviewed in this paper, as well as results from the search 
for millisecond pulsations in X-ray flux of the radio-quite 
neutron star RX J1856.5-3754. 

Keywords X-rays • Neutron stars • Millisecond pulsars 
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1 Introduction 

Millisecond pulsars (MSPs) significantly differ in prop- 
erties from other (ordinary) radio pulsars. First of all, 
MSPs possess very short and stable spin periods, P < 50 
ms, with extremely small period derivatives, P < 10~ 18 
s s _1 . These two main parameters separate MSPs from 
the majority of other pulsars, as illustrated in the P-P 
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Fig. 1 P-P diagram for about 1,500 radio pulsars (dots). 
Millisecond pulsars are located in the lower-left corner of the 
diagram. The pulsars from Table 1 and 3 are shown with red 
dots. Straight lines correspond to constant values of pulsar 
characteristic age r = 10 3 , 10 6 and 10 9 yr, surface magnetic 
field B SUTf = 10 8 , 10 10 , 10 12 and 10 14 G, and spin-down en- 
ergy E = 10 30 , 10 33 and 10 36 erg s" 1 . 



diagram 1 shown in Figure 1. According to the conven- 
tional pulsar magnetic-braking model, MSPs are very old 
neutron stars, with characteristic ages r = P/2P ~ 0.1 
10 Gyr, and low surface magnetic fields B sur { — 3.2 x 
10 19 (PP) 1 / 2 < 10 10 G (see Fig. 1). They were presum- 
ably spun up by angular momentum transfer during a 
mass accretion phase in binary systems 2 (e.g., Alpar et 
al. 1982), and their low values of i? sur f could be explained 

1 Based on the pulsar catalog provided by the Australia 
Telescope National Facility (Manchester et al. 2005) and 
available at http://www.atnf.csiro.au/research/pulsar 

2 It is why MSPs are often called "recycled" pulsars. 
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Table 1 Main parameters of eight MSPs 



PSR 


P 


d 


T 


logE 




(ms) 


(kpc) 


(Gyr) 


(ergs-) 


B1937+21 


1.56 


3.57 


0.24 


36.0 


B1957+20 


1.61 


2.49 


2.24 


35.0 


J0218+4232 


2.32 


2.67 


0.48 


35.4 


B1821-24 


3.05 


3.09 


0.03 


36.3 


J0030+0451 


4.87 


0.32 


7.71 


33.5 


J2124-3358 


4.93 


0.27 


6.01 


33.6 


J1024-0719 


5.16 


0.39 


> 27.25 


< 32.9 


J0437-4715 


5.76 


0.14 


6.51 


33.5 



by the Ohmic and/or accretion- induced decay of mag- 
netic field (see Cumming 2005 and references therein). 
Discoveries of several accretion-driven binary MSPs, first 
of all the famous pulsar SAX J1808.4-3658 with P ~ 2.5 
ms, support this hypothesis (see, e.g., Wijnands 2004 for 
a review on accreting MSPs). 

Since the discovery of the first fast rotating pulsar 
B1937+21 by Backer et al. (1982), MSPs have been ex- 
tensively searched for and studied in radio domain. Cur- 
rently, about 130 MSPs are known (Manchester et al. 
2005). Outside the radio band, as MSPs are intrinsically 
faint at optical wavelengths and most of them (<~ 80%) 
reside in binary systems with optically brighter white 
dwarf companions, the soft X-ray energy range (~ 0.1-10 
keV) is the main source of information on these pulsars. 
The detection of pulsed X-ray emission from the bright- 
est (and nearest) MSP J0437-4715 with ROSAT (Becker 
& Triimper 1993) initiated a series of dedicated X-ray 
observations of these intriguing objects in 90's with this 
satellite and also ASCA, BeppoSAX and RXTE. Later 
on this observational "relay" has been continued with 
Chandra and XMM-Newton. So far, firm X-ray detec- 
tions have been reported for about three dozens of iso- 
lated (solitary, or non-accreting if in binaries) MSPs. The 
majority of these objects are located in globular clus- 
ters (mostly 47 Tuc — see Bogdanov et al. 2006). This 
paper mainly concentrates on eight MSPs, listed in Ta- 
ble 1, with available detailed information on properties 
of detected X-ray emission. Besides P and t, Table 1 
gives estimates on the distances to these objects, d (in- 
ferred from either pulsar parallaxes or dispersion mea- 
sures and the NE2001 galactic electron density model 
by Cordes & Lazio 2003) 3 , and on their spin-down en- 
ergies, E = 4ir 2 IP~ 3 P (assuming a standard neutron 
star moment of inertia, I = 10 45 g cm 2 ). Note that for 
PSR J1024-0719 the intrinsic period derivative is not 
well determined, P < 3 x 10~ 21 s s -1 (Hotan, Bails & 



3 For the MSPs of the second group in Table 1 (with d < 1 
kpc), the distance estimates obtained from the pulsar paral- 
laxes and those inferred from their dispersion measures agree 
well with each other (see Hotan, Bails & Ord 2006 and Lom- 
men et al. 2006). 



Ord 2006), that results in the lower and upper limits on 
r and E, respectively. 

Generally, X-ray emission from radio pulsars consists 
of two different components, thermal and nonthermal, 
generated on the neutron star surface or in its vicin- 
ity. The nonthermal component is usually described by 
a power-law (PL) spectral model and attributed to ra- 
diation produced by synchrotron and/or inverse Comp- 
ton processes in the pulsar magnetosphere, whereas the 
thermal emission can originate from either the whole sur- 
face of a cooling neutron star or small hot spots around 
the magnetic poles (polar caps; PCs) on the star sur- 
face, or both. As predicted by virtually all pulsar models, 
these PCs can be heated up to X-ray temperatures (~ 1 
MK) by relativistic particles generated in pulsar accel- 
eration zones. A conventional assumption about the PC 
radius is that it is close to the radius within which open 
magnetic field lines originate from the pulsar surface, 
~ [2tt R^/cP] 1 ' 2 ~ 2[P/5ms]- 1 / 2 km (for a neutron 
star radius i?NS = 10 km). In case of MSPs, the entire 
surface at a neutron star age of ~ 1 Gyr is too cold, < 0.1 
MK, to be detectable in X-rays (although it may be seen 
in the UV/FUV band — see Sec. 6). Therefore, only 
nonthermal (magnetospheric) and/or thermal PC radi- 
ation is expected to be observed in X-rays from these 
objects. In addition to the radiation produced by MSPs 
themselves, nonthermal emission from pulsar-wind nebu- 
lae (PWNe) associated with MSPs moving at supersonic 
velocities ( > 100 km s -1 ) through interstellar medium 
may be detected. Another source of nonthermal X-ray 
radiation generated in binary systems can be an intra- 
binary shock formed where the pulsar wind and mat- 
ter from the stellar component collide (Arons & Tavani 
1993), although such a component would be hardly sepa- 
rated from radiation produced by the pulsar itself (unless 
properties of the nonthermal emission varies with orbital 
phase) . 

The MSPs in Table 1 belong to two distinct groups: 
those which emit almost pure nonthermal radiation and 
those with a predominantly thermal PC component. The 
next two Sections provides details on the X-ray proper- 
ties of these objects. Section 4 is devoted to three X- 
ray emitting MSPs whose properties remain uncertain. 
X-ray PWNe powered by MSPs are briefly discussed in 
Section 5. A summary is given in Section 6. 



2 Nonthermally emitting MSPs 

Of the MSPs listed in Table 1, these are PSRs B1937+21, 
B1957+20, J0218+4232, and B1821-24. They are char- 
acterized by large values of spin-down energy, E > 10 
erg s _1 , and shorter spin periods, P < 3.1 ms. PSR 
J0218+4232 is also one of few 7-pulsars known (Kuiper 
et al. 2000). X-ray observations of these MSPs revealed 
their nonthermal spectra, albeit with different photon 
indices T. Figures 2-5 show the spectra detected with 
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Energy (keV) 

Fig. 2 X-ray spectrum of PSR B1821-24 as detected with 
the Chandra ACIS-S instrument (crosses) fitted with a PL 
model of F — 1.2 (dotted curves) plus a possible PC compo- 
nent (dashes). See also Becker et al. (2003). 

PSR B1937+21 {Chandra) 
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3 Same as in Fig. 2 for PSR B1937+21 and a PL model 
= 1.2 (see also Kuiper et al. 2006). 
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Fig. 4 Same as in Fig. 2 for PSR J0218+4232 (as detected 
with the XMM-Newton EPIC-MOS instruments) and a PL 
model of T = 1.1 (see also Webb, Olive & Barret 2004). 
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Fig. 5 Same as in Fig. 2 for PSR B1957+20 and a PL model 
of r = 1.9. (Note that the detected flux can contain a con- 
tribution from a possible intrabinary shock — see Sec. 2 and 
Stappers et al. 2003 for more details.) 



Chandra and/or XMM-Newton from these MSPs. Al- 
though the spectral data on these pulsars do not formally 
require any other component (in addition to a PL), one 
cannot exclude that there is also a contribution of ther- 
mal PC emission in the detected X-ray fluxes, as indi- 
cated in Figures 2-5. For PSRs B1821-24, J0218+4232, 
and B1937+21 the measured photon indices are r ~ 1.1- 
1.2. The spectrum of PSR B1957+20, was found to be 
much steeper, with r ~ 1.9 (although the X-ray radia- 
tion of this pulsar can consist of two component — see 
below). Note that these estimates on r are derived for 
the phase-integrated fluxes, whereas the observational 
data on PSRs J0218+4232 and B1937+21 indicate that 



the spectral slope may change with pulsar rotational 
phase (see Webb, Olive & Barret 2004 and Nicastro et 
al. 2004). However, much more sensitive observations are 
required to confirm this effect. Regarding the X-ray spec- 
trum of PSR B1821-24 (located in the globular cluster 
M28), Becker et al. (2003) speculated that there is a 
marginal evidence for an emission line around 3.3 keV 
(see also Fig. 2). If this feature is real 4 , it could be in- 
terpreted as cyclotron emission from an optically thin 

4 There is no indication of a feature at this photon energy 
in the spectra of PSRs B1957+20 and B1937+21 detected 
with the same Chandra ACIS-S instrument (Figs. 3 and 5). 
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Fig. 6 Pulsed profiles of three nonthermally emitting MSPs. 
Vertical dashed lines indicate phases of main radio pulses. See 
Takahashi et al. (2000), Rutledge et al. (2004), and Kuiper et 
al. (2002) for more details on timing analysis of X-rays from 
PSRs B1937+21, B1821-24, and J0218+4232 (respectively). 

corona above the pulsar provided its magnetic field is 
strongly different from a centered dipole. The estimated 
nonthermal (isotropic) luminosities 5 of these four objects 
in the 0.2—10 keV range, £ nonth , and the corresponding 
"nonthermal" efficiencies, r? nonth = L nonth /E, are given 
in Table 2. Note that a fraction of the X-ray emission 
detected from the eclipsing pulsar B1957+20, which is 
thought to ablate its companion in a close binary sys- 
tem with a 9.2-hr orbital period, may be due to an in- 
trabinary shock between the pulsar wind and that of the 
companion star. The Chandra data on this MSP showed 
an apparent (at a 99% confidence level) modulation of 
the X-ray emission detected at the pulsar's position with 
orbital phase (Stappers et al. 2003), with lowest and 
highest fluxes during and immediately after eclipse (re- 
spectively). Assuming that this flux modulation is gen- 
uine, Stappers et al. (2003) obtained a 50% estimate on 
the contribution of X-ray emission from the intrabinary 
shock in the total flux detected from the pulsar. This 
fraction is accounted for in deriving the nonthermal lu- 
minosity and efficiency for PSR B1957+20 given in Ta- 
ble 2. 

Another prominent feature of the nonthermal emis- 
sion from three of these MSPs is the shape of the X- 
ray pulsed profiles, with strong and narrow main pulses 
and large pulsed fractions / p ranging from about 65% for 

5 X-ray luminosities of the pulsars discussed in Sections 2 
and 3 are derived for the distances given in Table 1. 



PSR B0218-4232 up to nearly 100% for PSRs B1937+21 
and B1821-24. Figure 6 presents the pulsars' light curves. 
Nonthermal pulsed emission from these MSPs was also 
detected up to about 20 keV with RXTE (Rots et al. 
1998; Kuiper, Hermsen & Stappers 2004; Cusumano et 
al. 2003). The main X-ray and radio peaks of these pul- 
sars were found to be nearly aligned in phase (based on 
both Chandra and RXTE data) 6 . This suggests that 
nonthermal photons emitted from these MSPs in the ra- 
dio and X-ray bands are generated in the same (or close) 
zones, although it still remains unclear where those zones 
are located: closer to the neutron star surface (as sug- 
gested by the polar-cap model of Harding, Usov & Mus- 
limov 2005) or near the pulsar light cylinder (according 
to the outer-gap model — e.g., Cheng, Ho & Ruderman 
1986, Romani 1996). Pulsations of the X-ray flux emitted 
by PSR B1957+20 have to be detected yet 7 . 



3 Thermally emitting MSPs 

This is the second group of objects in Table 1 which 
consists of PSRs J0030+0451, J2124-3358, J1024-0719, 
and J0437-4715. These MSPs are less energetic in terms 
of E, have longer spin periods and are much closer to 
Earth. X-ray spectra of the two brightest members of 
this group, PSRs J0030+0451 and J0437-4715 cannot 
be fitted with a single PL model, whereas for the other 
two pulsars such a fit yields too large photon indices, 
r ~ 3-4, and the obtained estimates on hydrogen col- 
umn density, njj — (1-2) x 10 21 cm~ 2 , towards these 
objects greatly exceed those inferred from independent 
measurements (Zavlin 2006). A broken-PL model, sug- 
gested by Becker & Aschenbach (2002) to interpret the 
spectra of PSRs J0030+0451 and J0437-4715, results in 
unrealistic estimates on 7Jh- A model involving a PL and 
a simple one-temperature thermal component also faces 
the same problem. 

Analyzing ROSAT data on PSR J0437-4715, Zavlin 
& Pavlov (1998) suggested a thermal model implying 
a nonuniform temperature distribution over PCs. These 
authors discussed that relativistic particles bombarding 
magnetic poles of an MSP could heat a region larger 
than the conventional PC size (see Sec. 1) because the 
low magnetic field of the pulsar does not prevent the 
released heat from propagating along the neutron star 
surface. The applied model assumes two identical PCs 
on magnetic poles covered with a weakly magnetized hy- 
drogen atmosphere (Zavlin, Pavlov & Shibanov 1996). 
It also takes into account the GR effects (redshift and 
bending of photon trajectories near the neutron star sur- 
face). In this model the thermal radiation depends on 

6 Although Takahashi et al. (2001) found from ASCA data 
that the main X-ray and radio pulses of PSR B1937+21 are 
separated by A(p — 0.5 in phase. 

7 No results from the XMM-Newton observation of this 
MSP conducted in October 2004 have been reported yet. 
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PC temperature and radius, the neutron star mass-to- 
radius ratio, which determines the GR effects, and the 
star geometry (the viewing and magnetic angles, £ and 
a, respectively). The sketch shown in Figure 7 illustrates 
the neutron star geometry and the effect of light bend- 
ing (see also, e.g., Zavlin, Shibanov & Pavlov 1995 for 
more details on the GR effects on pulsar PC emission). 
The nonuniform temperature was approximated with a 
step-like function, referred as PC "core" and "rim" . Ob- 
servations of PSR J0437-4715 with Chandra (Zavlin et 
al. 2002) and XMM-Newton (Zavlin 2006) showed that 
such a model, supplemented with a PL component of 
r ~ 2.0, fits well the spectrum of PSR J0437-4715 up 
to 10 keV and yields reasonable pulsar parameters as 
well as riH- The thermal model provides the bulk of the 
pulsar's X-ray flux, whereas the PL component prevails 
only at photon energies E > 3 keV. A similar (nonuni- 
form PCs plus PL) model works also well on the spec- 
tral data of PSR J0030+0451. Figures 8 and 9 demon- 
strate this. The inferred parameters of the thermal com- 
ponents are: T™ re ~ 1.4 and 2.1 MK, ~ 0.5 and 
0.8 MK, R c p °J c ~ 0.4 and 0.1 km, i?™ n ~ 2.6 and 1.4 
km, for PSRs J0437-4715 and J0030+0451, respectively 
(in the latter case the PL index was fixed at r = 1.5 
because of a poorer data statistics at higher energies). 
For PSRs J2124-3358 and J1024-0719 the same model 
yields PC parameters close to those mentioned above, 
although low quality of the observational data on these 
two objects at E > 3 keV allows one to put only upper 
limits on intensity of the nonthermal component (Zavlin 
2006). The bolometric luminosities of one PC, and 
"PC efficiencies", r; pc = L^JE, can be found in Ta- 
ble 2, together with the estimates on L nonth and r) nonth 
for the nonthermal component (or corresponding la up- 
per limits). Upper limits on the luminosity of a possible 
thermal PC component in the X-ray fluxes of the four 
nonthermally emitting pulsars (Sec. 2) are also given in 
Table 2. Note that a more detailed study of possible 
thermal emission from these four pulsars is hampered 
by large distances (a few kpc) and, hence, strong inter- 
stellar absorption (uh ~ [2-4] x 10 21 cm -2 ) towards these 
objects. 

As alternative to the PC-plus-PL model of the X-ray 
emission from PSR J0437-4715, Bogdanov, Grindlay & 
Rybicki (2006) suggested a purely thermal interpretation 
in which the harder X-ray spectral tail (at E > 3 keV) is 
a result of the inverse Compton scattering of soft ther- 
mal PC photons by energetic electrons/positrons in an 
optically thin thermal layer of a temperature kT c ~ 150 
keV located presumably in the pulsar's magnetosphere. 
However, as existence of such a thermal layer in a pulsar 
magnetosphere does not seem to be physically supported, 
this Comptonization interpretation remains rather ques- 
tionable. 

X-ray emission from all these four thermally emit- 
ting MSPs is pulsed, with pulsed fraction ranging from 
about 35% up to 50% (Becker & Aschenbach 2002; Zavlin 



| ^ observer 




Fig. 7 Sketch illustrating bending of photon trajectories in 
a strong gravitational field near the surface of a neutron star 
with rotational and magnetic axes 17 and m, and viewing 
and magnetic angles £ and a (respectively). 
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Fig. 8 X-ray spectra of PSR J0437-4715 as detected with 
the XMM-Newton EPIC-pn and MOS instruments (crosses). 
The solid curves show a best fitting model, PL (dots) plus a 
two-temperature PC component (dashes), "core" and "rim" 
(see Sec. 3). 



2006). The pulsed profiles of PSRs J0437-4715, J2124- 
3358, and J1024-0719 are rather similar in shape, with 
single broad pulses, whereas the light curve of PSR J0030 
+0451 exhibits two pulses per period separated by A(j) ~ 
0.5 in phase. This indicates that the geometry of PSR 
J0030+0451 (the angles £ and a) is different from those 
of the three others. For example, in the framework of 
the conventional pulsar model with the magnetic dipole 
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Table 2 X-ray luminosities of eight MSPs 



PSR LogL nonth Log L£p, Logr? nonth Log?7 pc 

(ergs" 1 ) (ergs- ¥ ) 



B1937+21 32.8 < 31.5 -3.2 < -4.5 

B1957+20 31.5 < 30.4 -3.5 < -4.6 

J0218+4232 32.6 < 31.0 -2.8 < -4.4 

B1821-24 32.7 < 32.0 -3.6 < -4.3 



J0030+0451 29.8 30.3 -3.7 -3.2 

J2124-3358 < 29.0 30.0 < -4.6 -3.6 

J1024-0719 < 28.9 29.3 > -3.6 

J0437-4715 29.7 30.2 -3.8 -3.3 




Energy (keV) 

Fig. 9 Same as in Fig. 8 for PSR J0030+0451. 



at the neutron star center, PSR J0030+0451 can be a 
nearly orthogonal rotator (i.e., C — ot — 90°) with two 
pulses in its light curve being due to contributions from 
two PCs seen during the pulsar's rotation. For the other 
MSPs, the bulk of the detected X-ray fluxes is expected 
to come mostly from one PC. However, the X-ray pulses 
of PSRs J0437-4715, J1024-0719, and J2124-3358 are 
clearly asymmetric, with a longer rise and a faster de- 
cay for the former two MSPs and with the opposite be- 
havior for the latter pulsar. None of these shapes can 
be explained by a simple axisymmetric temperature dis- 
tribution. A feasible interpretation is that the observed 
asymmetry in the pulsed profiles of these three MSPs 
is caused by contribution of the nonthermal component 
whose peak is shifted in phase with respect to the pulse 
of the thermal emission. Results of the energy-resolved 
timing and phase-resolved spectroscopy on the XMM- 
Newton data on PSR J0437-4715 support this explana- 
tion (Zavlin 2006). For PSR J2124-3358, there may be 
an alternative interpretation of the shape of the pulsar's 
light curve: the steeper rise and longer trail could be 
caused by relativistic effects (in particular, the Dopplcr 
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Fig. 10 Pulsed profiles of thermally emitting MSPs ob- 
tained from XMM- Newton observations (Zavlin 2006). The 
vertical dashed line in the upper panel indicates phase of a 
radio pulse of PSR J0437-4715 (see See. 3 for details). 



boost) in fast rotating pulsars (Braje, Romani & Rauch 
2000), although it then should be understood why these 
effects arc not seen in the pulse profiles of the other three 
MSPs whose the spin periods close to that of PSR J2124- 
3358 (Table 1). Another way to attempt to explain the 
observed light curves is to involve a model with a decen- 
tcrcd magnetic dipole (or another magnetic field config- 
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uration — see the presentation by S. Zane at this con- 
ference), what has not been done yet for interpreting 
pulsed X-ray emission from MSPs. In any case, to pro- 
duce a pulsed fraction at a level of f p > 30%, thermal 
emission has to be anisotropic, as predicted by neutron 
star atmosphere models, to counteract the effect of light 
bending near the star surface on pulsations of the ther- 
mal flux (Zavlin, Shibanov & Pavlov 1995). So far, the 
difference in phases between radio and X-ray pulses was 
determined only for PSR J0437-4715 based on Chandra 
data (Zavlin et al. 2002; see also Fig. 10). As the bulk 
of the X-ray flux from this pulsar is of the thermal PC 
origin, the small difference in phases of the radio and 
X-ray peaks suggests that the pulsar's radio emission is 
generated close to the neutron star surface, unless the 
effects of field-line sweepback and/or aberration cancel 
the travel-time difference. 



4 More X-ray emitting MSPs 

Three pulsars, J0034-0534, J0751+1807, and J1012+5307, 
are examples of MSPs firmly detected in X-rays but 
whose radiative properties remain uncertain. X-ray ob- 
servations of these objects were not long enough to pro- 
vide data suitable for discriminating among various spec- 
tral models. Those data could be equally well fitted with 
a PL model of r ~ 1.7 or a thermal (blackbody) model of 
a temperature of ~ 3 MK (see Webb et al. 2004 and Za- 
vlin 2006). Estimates on the total flux are the only rather 
reliable X-ray characteristics available for these objects. 
Table 3 presents the main pulsar parameters 8 together 
with the estimated X-ray luminosities in the 0.2-10 keV 
range, Lx, and corresponding efficiencies, rjx = Lx/E. 
In terms of the spin-down power E, PSRs J0751+1807 
and J1012+5307 join the four pulsars exhibiting thermal 
PC emission (Sec. 3), while PSR J0034-0534 seems to be 
more energetic and closer to the first group of the non- 
thermally emitting MSPs (Sec. 2). On the other hand, 
PSR J0034-0534 is most "underluminous" , its luminos- 
ity, Lx — 0.5 x 10 30 erg s _1 , is as low as that estimated 
for PSR J1024-0719, the least energetic pulsar among 
the MSPs discussed in Sections 2-4. The X-ray efficiency 
of PSR J0034-0534, i) X ^ 1.6 x 10" 5 , is smaller at least 
by an order of magnitude than those measured for the 
other MSPs (Tables 2 and 3). However, it is worth men- 
tioning that the estimate on E given in Table 2 for PSR 
J0034-0534 assumes that the intrinsic period derivative 
of the pulsar is as directly measured in radio observa- 
tions, i.e., without accounting for the Shklovskii effect 
on P due to pulsar proper motion (Shklovskii 1970). The 
proper motion of PSR J0034-0534 has not been deter- 
mined yet. Still, it is reasonable to assume that the trans- 

8 The distance to PSR J0751+1807 is derived from the pul- 
sar's parallax (Nice et al. 2005), whereas the other two esti- 
mates are obtained from the pulsars' dispersion measures and 
the NE2001 model. 



verse component of the pulsar's velocity may be about 
100 (c?/0.54kpc) 1 / 2 km s _1 (what is not unusual for a 
binary MSP — see Lommen et al. 2006). Then, its in- 
trinsic P and, hence, E would be reduced by a factor of 
10, making the pulsar's X-ray efficiency similar to those 

erg s" 



found for the MSPs with E < 10 34 — 



5 MSPs and their nebulae 

Relativistic pulsar winds, which carry away pulsar ro- 
tational energy, interact with ambient medium, that is 
expected to form PWNe detectable at different wave- 
lengths. In X-rays, about 30 PWNe are currently known 
(see Kaspi, Roberts & Harding 2006, Gaensler & Slane 
2006 and Kargaltsev & Pavlov 2006 for reviews, as well 
as electronic PWN catalogs 9 ), thanks mainly to the 
superb spatial resolution of Chandra (see the presen- 
tation by Weisskopf et al. at this conference for a re- 
view). The observed X-ray PWNe are diverse in prop- 
erties. Some of them are of a torus-like structure with 
jets along the symmetry axis (which apparently coincides 
with the pulsar's spin axis, as suggested for the young 
Crab and Vela pulsars — Weisskopf et al. 2000, Pavlov 
et al. 2003). Several others have a cometary-like shape 
caused by the pulsar motion. The "Mouse" PWN pow- 
ered by PSR J1747-2958 (Gaensler et al. 2004) is con- 
fined within a bow-shaped boundary without a shell-like 
structure. X-ray PWNe associated with PSR B1757-24 
(Kaspi et al. 2001), J1509 -5859 and J1809-1917 (San- 
wal, Pavlov & Garmire 2005), and the Geminga pulsar 
(Caraveo et al. 2003; De Luca et al. 2006; Pavlov, Sanwal 

6 Zavlin 2006) have elongated structures that look like 
"trails" (or "wakes") behind the moving pulsars. These 
tails could be pulsar jets confined by toroidal magnetic 
fields, or they could be associated with shocked relativis- 
tic wind confined by the ram pressure of the surround- 
ing interstellar medium. Two X-ray nebulae produced by 
the MSPs, B1957+20 (Stappers et al. 2003) and J2124- 
3358 (Hui & Becker 2006), belong to the latter group of 
PWNe. These are shown in Figures 11 and 12. Besides, 
optical observations revealed H Q bow-shocks surround- 
ing these two pulsars (sketched with dashed curves in 
Figs. 11 and 12). Both X-ray PWNe have a form of a 
20"-long tail, although there is a striking difference be- 
tween their shapes. The H Q bow-shock and X-ray tail of 
PSR B1957+20 are fairly symmetric, with the symmetry 
axis being nearly aligned with the direction of the pul- 
sar's proper motion (Fig. 11). For PSR J2124-3358, both 
the bow-shock and the X-ray tail are bent and highly 
asymmetric with respect to the proper motion vector 
(Fig. 12). Gaensler, John & Stappers (2002) argued that 
such an unusual shape of a PWN could be caused by a 
combination of effects associated with anisotropy of the 

9 http: //www. astro .psu. edu/users/green/psrdatabase/psr cat .htm 
http : //www .physics .mcgill . ca/ ~pulsar/pwncat .html 
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Table 3 Three MSPs with estimates on X-ray flux 



PSR 


P 


d 


T 


log-E 


logLx 


log?7x 




(ms) 


(kpc) 


(Gyr) 


(ergs -1 ) 


(ergs -1 ) 




J0034-0534 


1.88 


0.54 


6.0 


34.5 


29.7 


-4.8 


J0751+1807 


3.48 


0.63 


7.7 


33.8 


30.3 


-3.5 


J1012+5307 


5.26 


0.41 


8.6 


33.4 


30.4 


-3.0 



pulsar wind and density nonuniformity of the surround- 
ing medium. Analyzing the properties of the X-ray P WN 
of PSR B1957+20, Stappers et al. (2003) speculated that 
the efficiency with which relativistic particles are acceler- 
ated in the pulsar's wind can be as high as those inferred 
for young pulsars with much stronger surface magnetic 
fields. Comparing ratios of X-ray luminosities of tail- like 
PWNe to pulsar spin-down energies may be considered 
as an additional evidence in favor of this conclusion: 
the ratio, Lx/E ~ 10~ 4 , estimated for the tail associ- 
ated with PSR B1957+20 is very close to that found for 
the nebula produced by the 16-kyr-old PSR B1757-24 
(Kaspi et al. 2001) and greatly exceeds those inferred 
for the Geminga's tail, ~ 4 x 10~ 6 (Pavlov, Sanwal & 
Zavlin 2006) and the Vela's southeast jet, ~ 1 x 10~ 6 
(Pavlov et al. 2003). Contrary to these two objects, the 
third MSP whose supersonic motion through the ambi- 
ent medium causes an H a bow-shock, J0437-4715, has 
not been found to power an X-ray nebula, suggesting a 
very low magnetic field in the expected PWN region (Za- 
vlin et al. 2002). No an X-ray PWN associated with the 
most energetic MSPs B1937+21 has been detected in a 
deep Chandra observation of this pulsar (Kuiper et al. 
2006). 



6 Summarizing remarks 



* 



f20°47'55" 

19 h 59 n 



59-36?9 
K.A. 



Fig. 11 Chandra image of PSR B1957+20 and its tail-like 
X-ray PWN. The dashed curve sketches the H a nebula asso- 
ciated with the pulsar. The arrow indicates the direction of 
the pulsar's proper motion. 



X-ray observations of the MSPs discussed in this paper 
and those located in the globular clusters (Bogdanov et 
al. 2006) revealed that most of them emit predominantly 
thermal PC emission with a luminosity of ~ 10 30 - 
10 31 erg s- 1 and "PC efficiency" of r/P c - lO^-lO" 3 . 
It is worth noting that these estimates on n pc are close 
to that derived for the old ordinary pulsar B0950+08 
(r ~ 17 Myr, P ~ 253 ms, E ~ 6 x 10 32 erg s -1 ), 
if c ~ 2.4 x 10~ 4 (Zavlin & Pavlov 2004), suggesting that 
heating mechanisms in ordinary and rapidly spinning 
pulsars may be quite similar. Comparing pulsar param- 
eters of the MSPs with detected X-ray radiation shows 
that the four MSPs with nonthermal radiation possess 
much large values of E, by at least a factor of 50, than 
those estimated for the thermally emitting pulsars. The 
numbers presented in Table 2 indicate that the larger 
is E, the higher is the "nonthermal" efficiency rj nonth , 
whereas the opposite tendency is apparent for rf c . A sim- 
ilar behavior of if c follows from results of the theoretical 



modeling of the PC heating by returning positrons pro- 
duced through curvature radiation and inverse Compton 
scattering (see Figs. 7 and 8 in Harding & Muslimov 
2001 and 2002, respectively, with account for the rela- 
tion E ~ T -["+ 1 ]/["- 1 ] = T ~ 2 , for the magneto-dipolc 
braking index n — 3). In addition, high luminosity of the 
nonthermal emission may be associated with large values 
of the magnetic field at the pulsar light cylinder (Saito et 
al. 1997), Bic = B surf [27ri?Ns/cP] 3 . Indeed, the nonther- 
mally emitting MSPs possess magnetic fields, B\ c ~ (0.3- 
1) x 10 6 G, close to that of the Crab pulsar (~ 1 x 10 6 
G) and exceeding those in the MSPs with thermal emis- 
sion (as well as in PSRs J0751+1807 and J1012+5307), 
B\ c ~ 0.03 x 10 6 G, at least by an order of magni- 
tude. This may indicate that the nonthermal radiation of 
MSPs is generated in emission zone(s) close to the light 
cylinder (as predicted by the outer-gap pulsar models) 
and that B\ c is an important parameter governing the 
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Fig. 12 Same as in Fig. 11 for PSR J2124-3358. 



magnetospheric activity. On the other hand, a model by 
Harding, Usov & Muslimov (2005) for acceleration and 
pair cascades on open field lines close to pulsar PCs de- 
scribes rather well the spectrum of PSR J0218+4232 ob- 
served from X-rays through 7-rays. Hence, at the present 
stage, neither of the two models (outer-gap or polar-cap) 
of the nonthermal emission of MSPs could be regarded 
as more favorable. Besides, one cannot rule out that the 
pulsar spin period also affects jy nonth — the MSPs of the 
first group have shorter periods than those of the other 
four pulsars (Table 1). Definitely, a bigger sample of X- 
ray emitting MSPs is required to draw more conclusive 
speculations. 

The results of Section 3 point out that PCs of MSPs 
are likely nonuniform, with temperatures decreasing by 
a factor of 2-3 from the PC center down to its edge. 
There have not been yet reliable calculations of the tem- 
perature distribution around the pulsar magnetic poles. 
Hence, sophisticated theoretical models of PC heating 
and temperature distribution are required for further in- 
vestigation of thermal X-ray emission from MSPs. 

As shown by Pavlov & Zavlin (1997), modeling of an 
X-ray pulsed profile can yield constraints on the pulsar 
mass-to-radius ratio as well as its geometry. To do that, 
more elaborated models of magnetospheric pulsed emis- 
sion and data of a better quality at higher energies are 
required to disentangle nonthermal and thermal compo- 
nents. 

Chandra and XMM-Newton provided remarkable de- 
tections of two double-neutron-star binary systems with 



MSPs as companions, J0737-3039 and J1537+1155 (Cam- 
pana, Possentti & Burgay 2004; McLaughlin et al. 2004; 
Pcllizzioni et al. 2004; Kargaltsev, Pavlov & Garmirc 
2006). Although the collected data are of a rather scanty 
statistics, one could conclude that the inferred X-ray 
spectra and luminosities of these two system are similar, 
suggesting the same mechanisms generating the detected 
radiation. Among various interpretations proposed, the 
most plausible one is that the observed emission con- 
sists of a thermal PC component emitted by the MSPs 
plus nonthermal X-rays from the interaction of the pul- 
sar winds with the neutron star companions. This in- 
terpretation explains a gap observed in the orbital de- 
pendence of the X-ray flux emitted by the highly eccen- 
tric J1537+1155 system (Kargaltsev, Pavlov & Garmire 
2006). To investigate these systems in more detail, a 
phase-resolved spectroscopy on X-ray data of a better 
quality is required to separate possible components in 
X-rays from these objects. Further observations are ex- 
pected to constrain both the properties of the pulsar 
winds and those of the pulsars themselves in these (and 
other) double-neutron-star systems. 

So far, only one MSP, J0437-4715, has been observed 
and detected in the UV/ FUV band (Kargaltsev, Pavlov 
& Romani 2004). The shape of the inferred spectrum 
suggests thermal emission from the whole neutron star 
surface of a surprisingly high temperature of about 0.1 
MK. A heating mechanism should be operating in a Gyr- 
old neutron star to keep its surface at such temperature. 
As discussed by Kargaltsev, Pavlov & Romani (2004), 
among several possible mechanisms (both internal and 
external), chemical heating (Reisenegger 1995) and fric- 
tional heating (e.g., Cheng et al. 1992) of the neutron 
star core are the most plausible options. To understand 
thermal evolution of neutron stars, more MSPs (e.g., the 
other three pulsars discussed in Sec. 3) and close ordi- 
nary old pulsars (PSR B0950+08 is one of the best can- 
didates — Zavlin & Pavlov 2004) should be observed in 
the UV/FUV band. The example of PSR J0437-4715 
shows that such an observational program would be fea- 
sible. 

As a final remark, it is worth mentioning that so far 
all known isolated (non-accreting) MSPs were first dis- 
covered in radio. In X-rays a first attempt to find mil- 
lisecond pulsations from an isolated compact object was 
done by Zavlin & Pavlov (2006) who suggested that the 
most famous (and intriguing) member of the group of the 
"dim" isolated, radio-quite, and thermally emitting neu- 
tron stars discovered with ROSAT (see the presentation 
by F. Haberl at this conference), RX J1856. 5-3754, could 
be an MSP (see Pavlov & Zavlin 2003 for a discussion). 
One of the reasons to look for millisecond pulsations was 
the fact that the X-ray flux detected from this object 
in deep observations with Chandra and XMM-Newton 
revealed no pulsations at periods above 20 ms (the pe- 
riod range to which those observational data were sensi- 
tive) , with a stringent upper limit on the pulsed fraction 
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/ p < 1.3%. To check the MSP hypothesis, a deep XMM- 
Newton observation of RX J1856. 5-3754 in an instru- 
mental model with a temporal resolution of 0.03 ms was 
conducted in April 2004. Despite a very large number of 
source counts collected in this observation (~ 2 x 10 5 ), 
no significant pulsations were found at periods in the 
range P = 1-20 ms. The derived la upper limit on the 
pulsed flux for this period range is / p < 2.1%. This non- 
detection virtually excludes the MSP interpretation for 
RX J1856. 5-3754 and supports alternative hypotheses 
of a very special neutron star's geometry and/or a very 
long spin period (P > 10 hr). 
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